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ABSTRACT 

In galaxies like the Milky Way, cold 10^ K) gas ejected from the disc by stellar 
activity (the so-called galactic-fountain gas) is expected to interact with the virial- 
temperature 10^ K) gas of the corona. The associated transfer of momentum be- 
tween cold and hot gas has important consequences for the dynamics of both gas 
phases. We quantify the effects of such an interaction using hydrodynamical simula- 
tions of cold clouds travelling through a hot medium at different relative velocities. 
Our main finding is that there is a velocity threshold between clouds and corona, of 
about 75kms~"^, below which the hot gas ceases to absorb momentum from the cold 
clouds. It follows that in a disc galaxy like the Milky Way a static corona would be 
rapidly accelerated: the corona is expected to rotate and to lag, in the inner regions, 
by ^ 80 — 120 kms~^ with respect to the cold disc. We also show how the existence of 
this velocity threshold can explain the observed kinematics of the cold extra planar 
gas. 

Key words: hydrodynamics turbulence ISM: kinematics and dynamics Galaxy: 
kinematics and dynamics Galaxy: structure intergalactic medium 



1 INTRODUCTION 

Current cosmological models suggest that galaxies are em- 
bedded in extended virial-temperature atmospheres called 
coronae. These coronae are thought to contain a significant 
fraction (more than 50%) of the baryons in the Universe 
on the basis of a combination of big-bang nucleosynthe- 
sis theory and observations of the cosmic microwave back- 
sround (e.^- tF ukugita et al. 1998; Fukugita & Peebles 2004; 
iKomatsu et al.ll2009l ). In clusters of galaxies the baryons as- 
sociated with the hot intra-cluster medium are dir ectly de- 
tecte d in X-rays through their free-free emission ([Sarazinl 
l20Q9l ). whilst in lower density environments, such as the Lo- 
cal Group, the hot medium must be characterised by a much 
lower X-ray surface brig htness, below the sensitivity of the 
current instrumentation ( Anderson Bregmanllioiol V 

The exist ence of coronal gas around our Galaxy was hy- 
pothesised bv lSpitzeil (| 19561 ) to provide the necessary pres- 
sure confinement to cold clouds of interstellar gas that lie 
far from the Galactic plane. With the advent of UV spec- 
troscopy, highly ionised material (in particular C iv,0 V and 
O vi) was detected along lines of sight to distant UV sources 
that pass close to high- and intermediate-velocity clouds 



(hvcs and ivcs, see ISembach et al ] l2003l : iFox et al.ll2006l ). 

These clouds are not massive enough to be confined by self- 
gravity and the detected species are believed to lie at the 
interface between the coronal gas with T > 10^ K, and the 
HVCs/lVCs, which are composed by cooler and partly neutral 
gas. 

In star-forming disc galaxies, high-sensitivity Hi ob- 
servations have revealed that roughly 10% of the neutral 
gas is outside t he thin dis c, at distances of a few kilopar- 
secs (fOosterloo et al.l l2007l : iFraternali 2009 , and references 
therein). The origin of this extra-planar gas is not com- 
pletely understood but it is now widely accepted that a 
large fraction of it is produced by supernova- powered bub- 
bles which drive neutral and partly ionised gas out of the 
galactic disc. This gas is expected to travel through the 
halo and eventually to fall back to the disc in a time- 
scale of ~ 80 — 100 M yr, a mec hanism which is known 
as galactic fountain (Sh apiro k. Fie ld 1976; Bregman 198qI: 
iHouck Bregmanlll99d ). Thus, in the regions within a few 
kiloparsecs from the galactic disc a large number of galactic- 
fountain "cold" clouds (ivcs in the Milky Way) should be 
travelling through the ubiquitous coronal gas and an inter- 
action between the two phases is inevitable, with important 
consequences for the kinematics of both. 
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Unfortunately, the kinematics of the coronal gas is vir- 
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tually unconstrained: as mentioned, most of the evidence 
for this gas is indirect, and from the available data it is not 
possible to extract reliable kinematic information. A poten- 
tially useful direct detection comes from the observation of 
O VII absorption lines (which probe materi al at temperatures 
around 10^ K) towards bright quasars fe.g. lRasmussen et al.l 
However due to the poor spectral resolution of the 
curr ent X-ray instruments, these observations are im practi- 
cal (|Bregman fc Llovd-Davie3l2QQ7l : lYao et al.ll2Q08l V 

On the other hand, the kinematics of the extra planar 
cold gas has b een studied e xtensively with Hi and op- 
tical line data (|Swaters et al.l il 997: Fraternali et al. 20oJ; 
iHeald et al]l2007l : lOosterloo et al .1120071 : lKamDhuisll2QQ8i ). It 
was found that the extra planar gas rotates with a speed 
that decreases steadily with increasing distance from the 
mid-plane of the galaxy. A number of models have been 
proposed to explain this vert ical velocity gradient (e.g. 
Barnabe et aP l2006l: Kaufma niTet al.l l2006l : iMelioh etahl 



20091 : IStruck Smithli2009.. : .Marinacci et al.l 2010b V In par- 
ticular, it has been shown that the rotational gradient is 
larger than what would be expected for a collection of 
galactic-fountain clouds that travel through the halo re gion 
on pure ballistic trajectories (|Fraternali &: Binnevll2006l ). so 
the natural next step is to consider how clouds' trajectories 
are affected by interactions wit h the amb ient co rona. 

iFraternah fc BinnevI (12008", hereafter IfbosI ) considered 
two types of cloud-corona interaction: drag and accretion. 
The drag between the H i clouds and a supposedly lagging 
coro na was propose d, among others bv lCollins eHD (l2002l ) 
and iBarnabe et ahl ([2006), to explain the negative vertical 
rotation gradient of the extra planar gas. However, IfBOSI 
showed that this mechanism is not viable because of the re- 
sulting rapid acceleration of the corona, which causes the 
drag to vanish in less than a cloud orbital time. Thus, if the 
drag is as efficient as one would expect from purely mechan- 
ical arguments, the corona should rotate at a speed compa- 
rable to tha t of the disc, with a consequent vanishing of the 
drag. iFBOSi then investigated the possibility that H i clouds 
accrete some of the ambient gas that lies in their paths. 
They showed that if the ambient gas has a relatively low 
angular momentum about the galaxy's spin axis, its accre- 
tion onto H I clouds could explain the observed dynamics of 

the fountain clouds. 

In a preliminary study, iMarinacci et al.l (2010a, here- 
after iMlQ ) presented hydrodynamical simulations of cold 
clouds travelling through a hot medium with properties sim- 
ilar to the gas of a galactic corona. They found that radiative 
cooling can be efficient in the turbulent wakes of the clouds, 
allowing coronal material to condense in these wakes. There- 
fore, even though the original cloud is ablated, there is a net 
transfer of mass fro m the corona to the whole body of cool 
gas, supporting the lFBOSi framework. 

In this paper we extend the study of the interaction be- 
tween the clouds and the corona presented by iMld in two 
ways: (i) we use an upgraded version of the hydro-code tak- 
ing into account the variation in the mean molecular weight 
of the gas at T :^ 10"^ K (mainly due to hydrogen recombi- 
nation) and, more importantly, the different metallicity of 
the disc and of the coronal gas; (ii) we study a much wider 
range of relative cloud corona velocities. Our major find- 
ing is that there is a velocity threshold between cloud and 
corona below which the corona does not absorb momentum 
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Table 1. Parameters of the simulations. The initial tempera- 
ture and metallicity of the corona are Tcor = 2 x 10^ K and 
Zcor = 0.1 Zq, respectively, while for the cloud T^i = 104 K 
and Zci = Zq . All models are evolved up to 60 Myr. The grid 
size is 6114x1024 with spatial resolution ~ 2 x 2 pc for all the 
simulations. In the case of vlOO an additional run, with the same 
spatial resolution and grid size 6144x2048, was performed. 



from the cloud. We argue that the corona must be spinning 
fast enough for most clouds to move through the corona at 
a speed smaller than this velocity threshold. This prediction 
will be tested by future soft X-ray observations. 

The paper is organised as follows: Section [2] explains 
which simulations were run, and the results of these sim- 
ulations are analysed in Section [3] in Section |4l we discuss 
the implications of our results for the kinematics of both the 
corona and the clouds, while Section |5] sums up. 



2 SET UP OF THE HYDRODYNAMICAL 
SIMULATIONS 

T he init ial conditions of the simulations are similar to those 
of ImIoI : we simulate a cloud at T = 10"^ K with radius 
i^ci = 100 pc that moves at a speed vo through a station- 
ary and homogeneous backgrou nd at T — 2 x 10® K, which 
represents the corona. While in IMIOI the relative speed be- 
tween the cloud and the corona was fixed at vq — 75 kms~^, 
here we explore different relative speeds in the range 50 < 
Vo < 200 km s"-^. If these relative speeds are interpreted as 
differences in rotational velocities, then two extreme possi- 
bilities are encompassed by the simulations: (i) the case of 
a non-rotating corona (vq — 200kms~^) and (ii) a corona 
which is rotating with a speed close to that of the disc 
{vo — 50kms~^). As on ballistic trajectories the galacto- 
centric radii of the clouds would vary by less than 30%, and 
the clouds would reach heights above the galactic plane of 
a few kiloparsecs (see iFraternali Binnev l2006l ). the ex- 
pected dens ity va riations of the coronal gas are not strong. 
Thus, as in iMid we simplify our treatment by assuming a 
homogeneous corona and neglecting the effect of the gravi- 
tational field of the galaxy. We also neglect the self gravity 
of the gas, because the cloud masses are much smaller than 
their Jeans mass. At the beginning of each simulation the 
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cloud is in pressure equilibrium with its surroundings and 
the evolution of the system is followed for 60 Myr (signif- 
icantly longer than the 25 Myr spanned by most runs in 
In our standard simulations (see Table [T]) the density 
of the background is n — 10 ~^ cm~^. For each value of vq, 
we ran two simulations: one (dissipative) in which the gas is 
allowed to cool radiatively, another (adiabatic) without ra- 
diative cooling. The more realistic cases are those in which 
the gas is allowed to cool, but the adiabatic runs enable us to 
distinguish the effect of cooling (strongly dependent on the 
density, temperature and metallicity of the gas) from non- 
dissipative effects like mixing and hydrodynamical instabili- 
ties. Two additional sets of dissipative simulations have been 
run with half (low-density simulations in Table [T]) or twice 
(high-density simulations in Table [T]) the coronal density 
to assess the influence of this latter on the results. In all 
initial conditions the cloud's metallicity Zc\ takes the solar 
val ue ZfT) and the coronal metallicity is Z coY = O.IZq (see, 
e.g. lSembach et al.ll20Q3l : IShufl et al.ll2009l ). The metallicity 
of the fluid is allowed to vary from the initial conditions: 
in th is respect the present simulations differ from those of 
in which the metallicity was constant. This upgrade 
is important since the metallicity of a cloud ejected from 
the disc and that of the corona can differ by one order of 
magnitude or more, and the cooling rate of the gas critically 
depends on its metal content. In any case, we stress that the 
main results of lMld are conflrmed by the present investiga- 
tion: for instance, the global behaviour of iMld simulation 
with n — 10"^ cm~^ and uniform metallicity Z — 0.3 Zq 
is similar to ours with the same initial coronal density and 
cloud speed, but Zc\ — Zq and Zcor = 0.1 Zq. 



The simulations are performed with ECHO, a high- 
order, shock-capturing Eulerian hydrodynamical code: 
a detailed description and tests of it can be found 
in iDel Zanna &l Bucciantinil (l2002l ) and iDiT Zanna et al.l 
(l2007l ). In the current version of the code (echo-F-F) not 
only the metallicity evolves, but also the fluid's molecu- 
lar weight fi is a function of temperature. The variation of 
Z and fi is accounted for in the implementation of radia- 
tive co oling, which is modell e d acc ording to the prescrip- 
tion of [Sutherland Dopital (|l993l ) . The calculations are 
performed on a two-dimensional Cartesian grid with open 
boundary conditions imposed to all the sides of the compu- 
tational domain. Therefore, one of the dimensions perpen- 
dicular to the cloud's velocity has been suppressed, and we 
are in effect simulating a flow around an inflnite cylindrical 
cloud. The cylinder is moving perpendicular to its long axis 
and has initially a circular cross-section of radius i^ci- From 
the simulations we obtain quantities per unit length of the 
cylinder and we relate these to the corresponding quantities 
for an initially spherical cloud of radius Rd by multiplying 
the cylindrical results by the length |-Rci within which the 
mass of the cylinder equals the mass of the spherical cloud. 
This correction is included in the values of the cloud mass 
Mci quoted in Table [U We note that, due to pressure equi- 
librium in the initial conditions, for the same gas number 
d ensity , temperature and cloud size Md is higher here than 
m because here we account for the fact that fi depends 
on temperature, while in lMld fi is flxed at the coronal value. 



3 RESULTS OF THE SIMULATIONS 

We now describe the results of the hydrodynamical simu- 
lations, with emphasis on the mass and momentum trans- 
fer between the cloud and the ambient gas. While the two 
phases are unambiguously separated in the initial condi- 
tions (the cloud gas is at T = lO'^ K and the coronal gas 
is at T = 2 X 10^ K), at later times, as a consequence of 
turbulent mixing and radiative cooling, we expect to have 
a signiflcant amount of gas at intermediate temperatures. 
Therefore, at all times in the simulations, we deflne "hot" 
the gas at T > 10^ K and "cold" the gas at T < 3 x lO'^ K, 
identifying the hot phase with the corona and the cold phase 
with the cloud. We verifled that our results do not depend 
signiflcantly on the speciflc choice of these temperature cuts. 



3.1 Interaction between the cloud and the 
ambient medium 

Qualitatively, the evolution of the cloud in the present simu- 
lations is similar to that observed in the simulations of lMld . 
Due to the ram pressure arising from the motion, the cloud 
feels a drag which causes it to decelerate while the back- 
ground gains momentum. The body of the cloud is squashed, 
and behind it a wide turbulent wake forms, in which the 
coronal and cloud gas mix efficiently and radiative cooling 
is particularly effective, with the consequence that the cloud 
mass tend to increase with time. The wake has a complex 
velocity structure, with some fluid elements moving faster 
and other slower than the body of the cloud. 

This behaviour is apparent from Figs. [T] and |2l show- 
ing, respectively, density and velocity snapshots, at four 
different times, of the dissipative standard simulation with 
vo — 100 km s"^ (right panels) and, for comparison, of the 
corresponding adiabatic simulation (left panels). The dia- 
grams in the bottom panels clearly show that the cloud 
is decelerated by drag: at lOOkms"^ ^ 0.1 kpcMyr"^ the 
cloud would travel, with no drag, a distance of about 6 kpc 
in 60 Myr, but neither in the presence nor in the absence 
of radiative cooling does the main body of the cloud reach 
that distance. The structure of the body of the cloud and 
of the wake is strongly influenced by the presence of radia- 
tive cooling: the wake is more elongated and less laterally 
extended (except perhaps in the 60 Myr panel) when the ra- 
diative cooling of the gas is permitted. The extent of the 
velocity structure of the wake is greater, both laterally and 
in elongation, in the dissipative than in the adiabatic case at 
late times (see Fig. |2j. Just behind the cloud there is a re- 
gion in which the velocities reach their maximum values and 
this region is more compact (less elongated) when cooling 
occurs. In front of the cloud the coronal gas is pushed by the 
cloud and therefore also in this region the gas has positive 
velocity. However, in contrast to what happens behind the 
main body of the cloud, this region is more extended in the 
absence of cooling. 

Comparing the flnal conflguration with the initial con- 
ditions, we flnd that the cloud has increased its mass (when 
cooling is allowed) and that the corona has gained momen- 
tum. Quantitatively, the amount of mass and momentum 
transfer depends critically on radiative cooling. This can 
be seen in Fig. [S] which shows the distributions of mass 
(upper panel) and momentum (lower panel) as functions of 
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Figure 1. Density snapshots of the adiabatic (left panels) and dissipative (right panels) standard simulations with vq = 100 km s ^ 
(vlOO in Table 0). The time at which the snapshot has been taken is indicated in each panel. The initial position of the cloud centre is 
X = 0. 
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Figure 2. Same as Fig. but plotting the velocity along the cloud's direction of motion {va 



temperature after 50 Myr for the standard simulation with 
'L'o = lOOkms" and, for comparison, for its adiabatic coun- 
terpart. In both panels the histograms peak around the tem- 
peratures of the coronal gas and the initial cloud tempera- 
ture (T = 2 X lO^K and T = lO'^K, respectively), and the 
distribution around these peaks is rather narrow. From the 
bottom panel of Fig. O it is clear that there is an exchange 



of momentum between the cloud and the corona, but the 
corona acquires less momentum when the gas is allowed to 
cool, as can be inferred from the heights of the peaks at 
T > 10^ K. This can be better understood if we recall how 
the mass of the gas is distributed among the various phases. 
At the beginning, most of the gas is at the cloud and the 
corona initial temperatures, and only a small fraction of the 
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Figure 3. Distributions of mass (upper panel) and momentum 
(lower panel) as functions of the gas temperature after 50 Myr 
of evolution of the the adiabatic (shaded) and dissipative (long- 
dashed line) standard simulations with = 100 km s~^; the black 
short-dashed line is the initial mass or momentum distribution. 
The vertical grey dashed lines represent the temperature cuts 
used to define cold and hot gas phases. 



total mass is at intermediate temperatures. The momentum 
of the cloud is mostly transferred to this intermediate tem- 
perature component as a result of the mixing. Because this 
gas is in the temperature range T ~ 1 — 5 x 10^ K, in which 
the cooling function reaches its maximum, it can cool very 
effectively. The consequence of this cooling is a mass trans- 
fer from the corona towards the cold gas. The momentum 
removed from the cold cloud is thus retained by the cooling 
gas and never transferred to the coronal gas. If the cooling 
is not present, the process of condensation of the mixed gas 
cannot occur and therefore the cloud continues to transfer 
its momentum to the mixed gas and thus to the corona. 



3.2 Motion of the cold gas 

In order to quantify the deceleration of the cold cloud, in the 
simulations we monitor the evolution of its centroid velocity, 
defined as the total momentum of the cold gas in the cloud's 
direction of motion x divided the total mass of the cold gas. 
We com pare the centroid velocity with the analytic estimate 

(IfbosI ) 



where the characteristic drag time is given in terms of the 
cloud's initial mass Md, geometrical cross section a, and 
coronal density ph, by 



Mcl 



(2) 



In all our dissipative standard simulations, independent of 
the initial speed vo, after 60 Myr the centroid velocity is 
~ 0.75 Vo, significantly lower than the values predicted by 
equation ([T]), which treats the cloud as a rigid body. This dis- 
crepancy is due to a combination of cloud deformation, hy- 
drodynamical instabilities and condensation of cooling gas. 
This can be inferred from Fig. lU which shows the cen- 
troid velocity for the dissipative (stars) and adiabatic (plus 
signs) standard simulations with -uo = 200kms~^ (top), 
Vo — lOOkms"^ (centre) and vo — 75kms^^ (bottom). It is 
clear that the analytic formula (solid line in each diagram) 
is a good description of the general behaviour of cold gas 
motion over ~ 60 Myr only when the cooling is not allowed 
and the initial relative speed is sufficiently low (see bottom 
panel in the figure). When these conditions are not satisfied, 
the analytic formula fails after some time, mainly because 
Kelvin-Helmholtz instability strips some cold gas and the 
ram pressure progressively squashes the cloud increasing its 
cross-section. The net effect is an enhancement of drag ef- 
fectiveness (see equation |2|. 

Let us now compare the evolution of the centroid ve- 
locity in the presence and in the absence of cooling. First, 
it must be noticed that in normalised units {vo for veloci- 
ties and Rc\/vo for times, shown in the right and upper axes 
in the diagrams in Fig. |4|) the evolution of the centroid ve- 
locity is basically the same in all the adiabatic simulations 
(in the overlapping normalised time intervals). In particular, 
the deviation from the analytic estimate becomes apparent 
at roughly the same normalised time and it can be ascribed 
mainly to the flattening of the main body of the cloud. On 
the other hand, in the presence of cooling the deviation of 
the cloud speed from the analytic estimate occurs at roughly 
the same physical time (~ 20 Myr) for all explored initial 
speeds. In these dissipative cases the amount of gas that 
transfers from the hot to the cold phase has an appreciable 
impact on the centroid velocity, so we are no longer merely 
looking at the slowing of the original cloud, and the deceler- 
ation is dominated by condensation of cooling material that 
was originally at rest. 



3.3 Momentum of the hot gas 

We now turn to the momentum transferred to the hot gas. In 
the simulations, the total momentum of the hot gas is com- 
puted by summing, over the whole computational domain, 
the momentum of the cells containing gas at T > 10^ K. 
Figure |5] shows the amount of momentum, along the cloud's 
direction of motion x, acquired by the hot g as as a func- 
tion of time for the dissipative (stars) and adiabatic (plus 
signs) standard simulations with vo — 200, 100, 75kms~^. 
From these diagrams it is apparent that the hot gas gains 
less momentum in the presence than in the absence of ra- 
diative cooling. The ratio between the flnal coronal momen- 
tum in the dissipative and in the adiabatic simulations is 
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Figure 4. Evolution of the centroid velocity, along the cloud's 
direction of motion, of cold gas in the dissipative (stars) and adi- 
abatic (plus signs) standard simulations with vq = 200kms~^ 
(top panel), vq = lOOkms"^ (central panel) and vq = 75kms~^ 
(bottom panel). The solid lines are the theoretical predictions in 
the case of a quadratic drag (equation 0). 



~ 0.8 when vq — 200kms~^ (top panel), ~ 0.5 when vo = 
lOOkms"^ (middle panel), and ^0.1 when vq — 75kms~^ 
(bottom panel). In particular, when vq — 75kms~^ the mo- 
mentum transferred from the cold gas to the corona, after 
an initial phase of growth lasting for ~ 20 Myr, settles to a 
nearly constant value: from that moment onwards, the mo- 
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Figure 5. Evolution of momentum acquired by the corona as 
a function of time in the dissipative (stars) and adiabatic (plus 
signs) standard simulations with vq = 200, 100, 75 kms~^ (from 
top to bottom). For comparison, we note that the total initial 
momentum of the cloud is 2.4 x 10® (-uo/lOO km s~^). 



mentum transferred to the corona is consistent with zero. 
Thus, the simulations indicate the existence of a velocity 
threshold Vth below which the transfer of momentum from 
the cold clouds to the corona is suppressed. Given that the 
cloud keeps losing momentum (Fig. [4)) and the total momen- 
tum is conserved, it follows that the momentum lost by the 
cloud is retained by the coronal gas, which is progressively 
cooling onto the cloud's wake. 

The measure of the momentum transferred from the 



cold to the hot gas can be affected by the flow of coronal 
gas (and associated momentum) through the open bound- 
aries of the computational domain. This effect becomes non- 
negligible only in the last ~ 10 Myr of the simulations, so 
we take ~ 50 Myr as flnal time, as far as the measure of the 
coronal momentum is concerned (see also Fig.|3j. As an ad- 
ditional test, we also ran a simulation with the same initial 
condition as the standard simulation with vq = 100 kms~^, 
but with twice the number of grid points in the direction 
perpendicular to the cloud's direction of motion (keeping 
the resolution flxed and thus effectively doubling the size of 
the domain), flnding relative differences in the actual value 
of the coronal momentum within 8%. 

Since the density of the corona strongly influences the 
cooling rate of the gas, we tested whether the velocity 
threshold depends on this parameter with two additional 
sets of simulations with half and twice the above coronal 
density. In particular we ran two simulations with higher 
density (vo — 75, 100 km s~^) and three simulations with 
lower density (vo — 50, 75, 100 kms~^; see Table [T|. In a 
denser corona, due to the increased cooling rate, the value 
of the velocity threshold rises to 85kms~^. In the case of 
the less dense corona, we found that momentum is trans- 
ferred signiflcantly also when vo — 75kms~^, but when 
'L'o = 50kms~^ the momentum transfer is consistent with 
zero, so Vth ~ 50kms~^. This means that the actual value 
of the velocity threshold Vth is sensitive to the coronal den- 
sity, and therefore uncertain. However, provided that in 
Milky Way like galaxies the actual average coronal density 
a few kiloparsecs above the plane is in the explored interval 
0.5 — 2 X 10~^ cm~^, we can consider Vth ~ 50 — 85kms^^ 
as a fiducial range for the velocity threshold. 



4 IMPLICATIONS FOR GALACTIC 

CORONAE AND COLD EXTRA-PLANAR 
GAS 

4.1 Rotating coronae 

IfbosI showed that if a significant part of the momentum lost 
by H I clouds through ram-pressure drag is absorbed by the 
corona, the spin of the latter will increase on a time-scale 
that is short even compared to the orbital time (~ 100 Myr) 
of an individual cloud. If the corona did not rotate, the 
velocity differences between it and individual clouds would 
be dominated by the cloud's azimuthal motion, and being 
>200kms~^ would exceed the velocity threshold vth- Con- 
sequently, the corona would gain a significant fraction of 
the momentum lost by clouds, acquiring substantial angular 
momentum on a short time scale. Hence the rotation rate 
of the corona must be determined by the condition that 
the relative velocity of a typical cloud and the corona is 
<vth- Since clouds are tr avelling with an average ve rtical 
velocity (vz) ~ 40kms"-^ (|Marasco fc Fratern'aiill201lh . the 
rotational lag required for this condition to be satisfied is 



-uiag ~ V'^th ~ i'^z)'^- For the estimate of vth ~ 75kms~^ 
given above, we obtain -uiag ~ 65kms^^. Thus, we expect 
coronae to be quite rapidly rotating. 

To compute the difference in rotation velocity between 
the corona and the disc at a given height, one must subtract 
the value for -uiag determined above to the value of the rota- 
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Figure 6. Meridional density contours of a rotating isothermal 
(T = 2 X 10^ K) corona in the Milky Way potential for two values 
of the asymptotic rotational velocity (see equation The 
density is normalised to the central {R = z = 0) value and the 
values on the contours are 0.45, 0.4, 0.35, 0.3, 0.25, 0.2, 0.15 (the 
lowest contour is not visible in the bottom panel). 



tion velocity of the H I at that height. In the Milky Way, as- 
suming a vertical gradient fo r the H I of —15 kms~^ kpc~^ 
(IMarasco &: Fraternalill20Tll ). we find that this difference is 
~ 95kms~"^ at 2; ~ 2 kpc, corresponding to a coronal rota- 
tion velocity of ~ 125kms~^ for vq — 220kms~^. 

The fact that the corona is rotating, even at speeds 
lower than that of the galactic disc, has important conse- 
quences for its morphology. In general, the azimuthal ve- 
locity of the corona can depend on both distance from 
the centre R and height above the plane z, as is t he case of 
baroclinic distributions (Poincare Wavre theorem: lLebovit3 
^67; Tassoul 1978; Barnabe et al. 2006). The z dependence 
might be expected in analogy with the observed kinematics 
of the cold extra-planar gas, and also as a consequence of the 
momentum transfer between fountain clouds and hot gas. By 
recalling the procedure outlined above, the lag of the corona 
with respect to the disc should vary from ~ 70 to 120 kms~^ 
going from z — 0.5 kpc to 2; = 4 kpc, which is the typical re- 
gion where fountain clouds are located. This corresponds to 
a rotational velocity of ~ 150 — lOOkms"^, assuming again 
Vq — 220kms~^. However, given the large uncertainty on 
the properties of the corona, here we consider much simpler 
(barotropic) models of coronae, in which v^p — v^p{R). In 
particular, we construct isothermal (T = 2 x 10^ K) models 
of coronae in equilibrium in the axisy mmetric gravitational 
poten tial of the Milky Way, as given bv lBinnev Tremaind 
(|2008l ) . with rotation law 



v^{R) 



R 



(3) 



The rotation is null at the centre and raises up to the ter 
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minal value Vs'. the rise of the rotation velocity is steeper 
for smaller values of Rs. To match the expected rotational 
velocity range, in Fig.[6]the isodensity contours of two mod- 
els with Vs — lOOkms"^ (upper panel) and Vs — 150 km s^^ 
(lower panel) are presented. In both cases Rs — 1 kpc. The 
resulting gas distributions are peanut-shaped, with axis ra- 
tio of the isodensity surfaces in the range 0.5—0.7. Therefore, 
if galactic coronae rotate with speeds of the order of those 
predicted by our model, they will appear significantly fiat- 
tened when detected in future X-ray observations. We note 
that this peanut-shape d structure is a lso predicted in the 
general baroclinic case (|Barnabe et al.l l2Q06V 

The kinematic properties of the galactic coronae are im- 
portant also for their evolution. For instance, the thermal- 
stability properties of a corona can be significantly in- 
fiuen ced by the d istribution of specific angular momen- 
tum (lNiDotill2010l). with implications for the origin of the 
HVCs (|Binnev et al.ll2009l V However, the HVCs are typically 
thought to be relatively far from the disc, in a region not 
strongly affected by disc activities, while in the present work 
we have focused on the kinematics of the corona within a 
few kiloparsecs from the disc, where the galactic fountain 
is at work. It is not clear to what degree the kinematics of 
the hot gas far from the disc (which is likely related to the 
formation history of the galaxy) can be influenced by the 
kinematics of coronal gas near the disc. 



4.2 Vertical velocity gradient of the cold 
extra— planar gas 

The simulations presented here clearly indicate that clouds 
of cold gas ejected above the disc plane can lose angular mo- 
mentum via interaction with the hot medium. This mech- 
anism appears as a natural explanation for the observed 
vertical gradient in the rotation velocity of the extra planar 
gas of disc galaxies. Unfortunately, a comparison between 
the results of the simulations and the observations is not 
straightforward, mainly because of the idealised nature of 
the simulations, in which we do not include the galaxy's 
gravitational fleld and thus we do not follow the cloud's orbit 
above the plane. However, our results can be used to obtain 
a rough estimate of the expected vertical rotational gradient 
of the extra planar gas, using the following argument. Let 
us simplify the treatment by assuming that a cloud, ejected 
with a velocity orthogonal to the disc plane at some radius 
i?, falls back onto the plane, after a time of the order of 
100 Myr, at about the same radius R. Assuming a mean ve- 
locity for the ascendirig par t of the orbit (vz) — 40kms~^ 
([Marasco Fraternal] l2Qlll ) and a mean maximum height 
above the disc for the cold clouds {zmax) — 2 kpc, the time 
to reach the maximum height is 



^orb 



50 Myr. 



(4) 



Thus, neglecting for the moment any change in R along the 
orbit, the rotational- velocity gradient of the cold gas can be 
estimated as 



Az 



Av 



('2'max) 



(5) 



where Av — v(torh) — vo is the difference between the cen- 
troid velocity of the cold gas after torh and its initial velocity. 



In the hypothetical case of a static corona {vo = 200 km s"^) 
Av ^ 40kms~^ (see Fig. U, so 



Av^ 
Az 



-20 km s ^ kpc ^ , 



(6) 



while in the case of a spinning corona (such that vq — 
75 kms-M 



Avcp 
~A^ 



-7kms ^ kpc ^ 



(7) 



based on the results of the dissipative simulation v75. We 
note that the effect of radiative cooling is crucial for the 
estimate of the gradient: using the adiabatic vq — 75kms~^ 
simulation we would get instead 



Av^p 
~A^ 
In 



-3kms ^ kpc ^ 



(8) 



the Milky W ay the observed gradie nt is 
^ — 15 kms~^ kpc~^ ( Marasco Fraternalil l201lh . and 
similar values have been determined in external nearby 
galaxies. As mentioned, the above estimates neglect the 
variation of R along the orbit, which by conservation of 
angular momentum induces a variation in v^p. IfbosI showed 
that this effect alone accounts for ~ 50% of the rotational 
gradient of the extra-planar gas. Thus, by adding this 
latter contribution to the estimate in equation ([7|), we get 
remarkably close to the observed gradients. 

On the basis of the above discussion we can conclude 
that: (i) larger velocity gradients are expected with increas- 
ing differences in rotational velocity between the clouds and 
the corona, (ii) as the relative velocity decreases, the drag 
becomes an inefficient process to decelerate the clouds, (iii) 
the radiative cooling of the gas has the effect of increasing 
the gradient to values which are very close to those required 
by the observations, and therefore (iv) the observed gradi- 
ent can be explained only by a combination of cooling and 
drag. Over all, our results are consistent with the scenario de- 
veloped in 'FB08', in which the accretion of gas by fountain 
clouds is the cause of the lagging kinematics of extra planar 
gas in disc galaxies. 



5 SUMMARY AND CONCLUSIONS 

In this paper we studied the process of the momentum 
transfer between cold clouds ejected from the disc and the 
corona of disc galaxies like the Milky Way, which has a sig- 
niflcant impact on the kinematics of both the coronal and 
the cold extra planar gas. In particular, we presented two- 
dimensional hydrodynamical simulations of cloud-corona in- 
teraction. In the more realistic simulations the gas is allowed 
to cool radiatively, but, for comparison, we also considered 
the corresponding adiabatic cases. The simulations are simi- 
lar, but signiflcantly improved, with respect to those carried 
out by IMIQ : here we explored a range of relative speeds 
between the cloud and the hot medium (50 — 200kms~^), 
we followed the evolution of the system for longer time 
(60 Myr), we treated the metallicity as a dynamical variable 
(assigning initially to the cloud a metallicity ten times higher 
than that of the ambient gas) and we accounted for the vari- 
ation with temperature of the fluid's molecular weight /i. 
From the simulations we draw the following conclusions. 
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(i) There is a relative velocity threshold of about 
75kms~^ between the cloud and the ambient medium, be- 
low which the corona stops absorbing momentum. 

(ii) Radiative cooling is crucial for the existence of such 
velocity threshold: in adiabatic simulations the momentum 
transfer from the cloud to the hot gas is quite efficient at all 
explored relative speeds. 

(iii) The actual value of the velocity threshold for Milky 
Way like galaxies depends on the coronal density: less dense 
coronae have a lower velocity threshold. In particular, in 
the interval of coronal densities explored by our simulations 
vth ~ 50 — 85kms~^. 

The process studied here has an important consequence 
on the kinematics of the corona. The corona cannot be static, 
otherwise the momentum transfer from the cold clouds 
would be very effective, but it cannot be rotating at a speed 
close to the rotation speed of the disc, because at low rela- 
tive disc-corona (and thus cloud-corona) velocity the cooling 
hampers any momentum transfer. In a galaxy like the Milky 
Way we expect the corona to lag, in the inner regions, by 
~ 100 km s~^ with respect to the cold disc at z c::^ 2 kpc. Such 
a spinning corona must be characterised by a significantly 
flattened density distribution: based on simple isothermal 
models, we estimate values of the axis ratios of its isoden- 
sity surfaces in the range 0.5 — 0.7. 

Due to interaction with the corona, fountain clouds gain 
mass and decelerate: based on our calculations, this effect 
can account for — 7kms~^kpc~^ of the vertical gradient 
in rotational speed of the extra-planar gas in disc galaxies 
like the Milky Way. Such a value can be sufficient to reconcile 
purely ballistic fountain models with the observed kinemat- 
ics of the Hi extra planar gas, giving further support to the 
"fountain plus accretion" model of FBO& . 

In conclusion, we think that the presented results con- 
tribute quantitatively to refine an emerging consistent sce- 
nario in which galactic fountains, cold extra planar gas and 
galactic coronae are strictly interlaced. This scenario makes 
specific predictions that can be tested in the future with 
soft X-ray observations detecting the elusive galactic coro- 
nae and with improved H I measures of the cold extra planar 
gas. From the theoretical point of view, the next step will 
be the extension to three dimensions of the two-dimensional 
calculations here presented, and ideally one would like to 
realise "global" simulations, in which the entire galaxy is 
modelled, with supernova- driven fountain clouds orbiting 
through the hot coronal medium under the effect of the 
galactic gravitational field. 
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